Impact of defect layers on superconductivity in Bi-O-S compounds 
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Newly discovered Bi-S-0 compounds remain an enigma in attempts to understand their electronic 
properties. Recent study of Bi404S3 has shown it to be a mixture of two phases, Bi20S2 and 
3130283, the latter being superconducting. Using density functional theory, we explore the electronic 
structure of both the phases and the effect of the introduction of stacking faults. Our results 
demonstrate the S2 layers dope the BiS2 bands, these electrons are accumulated in the BiS planes. 
We show the introduction of defects in the stacking diminishes the doping which would suppress 
superconductivity. 



The recent discovery of superconductivity in two new 
compounds, 6140483 [1-4] and La02:Fi_a;BiS2 [5-8], has 
raised great interest. In fact it is a whole new class of 
superconductors that has been unveiled, as a number 
of BiS2-based compounds exhibit superconductivity. All 
those compounds share the same structure : a stacking 
of alternating BiS2 and spacer layers. Superconductivity 
is thought to arise from doping of the BiS2 layers. Both 
hole, Sri_2;La2;FBiS2 [ ], and electron, Lai_j,M2;OBiS2 
(M = Ti, Zr, Hf, Th) [.:■], doping result in supercon- 
ductivity. The versatility of the La(0,F)BiS2 system has 
been demonstrated by replacement of La with a range of 
lanthanide Ln^~^ ions [11-15] — a maximum Tc of 10.6 K 
is reported in La02,Fi_a;BiS2 for x = 0.5. However given 
the absence of a systematic study and the difficulties with 
sample preparation, it may still be possible that improve- 
ment in condition and control of disorder might result in 
achieving higher temperature superconductivity. 

Electronic structure calculations have mainly focused 
on the La02:Fi_a;BiS2 (x = 0, 0.5) materials partially due 
to uncertainties in the composition of the Bi-O-S super- 
conducting phase. These calculations indicated that the 
superconducting electrons are a mixture of Bismuth 6p 
and Sulphur 3p states [16, 17]. Electronic and thermo- 
dynamic properties of this group of materials are indeed 
exciting and enigmatic and suggestions range from spin 
fluctuations mediated superconductivity [IS] to proxim- 
ity to ferroelectricity and Charge Density Waves (CDW) 
instability [1')]. 

Despite the excitement, doubts still remain concern- 
ing the structure of the first superconductor, 6140483. 
However Phelan et al. just published an extensive study 
of the chemistry of these compounds and concluded that 
6140483 was actually a two-phased material [2U]. The 
two phases are Bi2082 and 6130283, the later being likely 
the superconducting one. 
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Composed of 6182, S2 and 6i202 layers, the latter be- 
ing a spacer layer (Figure 1), 6130283 is made of alter- 
nating 6182 and 82 layers, all separated by spacer layers, 
whereas in 612 082 only the two 6182 layers are present. 
Initial experimental results suggest that superconductiv- 
ity is suppressed by the inclusion of additional 6182 layers 
disrupting the alternate of 82 and 6182 sheets. Stacking 
faults play a significant role in superconducting proper- 
ties of several families of materials [-1]. For example, 
changes to the ground state have been observed in the 
presence of stacking faults in other systems, including 
modifications to superconductivity [22, 23]. 

In this paper we explore the electronic structure of the 
superconducting and non-superconducting phases of the 
6i-0-8 compounds. The role of introducing additional 
6182 layers in the superconductor 6130283 is also inves- 
tigated. Like in the La(0,F)6i82 systems we find that the 
6182 layers play a critical role in superconductivity. Our 
analysis indicates that the 82 sheets donate significant 
electron density to the 6182 bands. Those electrons are 
localised within 6i8 planes rather that through the 6182 
sheets. Interruption of the 6182, 82 stacking sequence in 
defect structures is shown to reduce electron donation to 
the 618 planes, suppressing superconductivity. 

6and structures and density of states for 6i2 0S2 and 
6130283 were calculated, using the SIESTA code [24, 25], 
implementing the generalized gradient approximation in 
the shape of the Perdew, 6urke, and Ernzerhof func- 
tional. It uses non-conserving pseudopotentials to re- 
place the core electrons, while the valence electrons are 
described using atomic-like orbitals as basis states at the 
double zeta polarized level. In the case of the two known 
structures, experimentally determined structural param- 
eters were used in our calculations [20] [28]. Spin polar- 
ized calculations were found to be unnecessary. 

In order to simulate how stacking faults perturb the 
electronic structure, we investigated materials where the 
relative number of 6182 and 82 layers was changed. In 
6130283, the stacking is 2 x 6182 / 61202 / 82 / 61202 
/ 2 x 6182 / Bi202 / 82 / 61202 (Figure 1). The se- 
quence is repeated twice because of symmetry : 82 layers 
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FIG. 1: Crystal structure of Bi20S2, Bi302S3, 6130283 with 
an added stacking fault, and two BI3O2S3 on top of each other 
with an added stacking fault, represented In Vesta [2()]. They 
are formed of a stacking of BIS2, S2 and BI2O2 layers, the 
latter being a spacer layer. These structures Implement four 
different 82/6182 ratios : 0, 1, 2/3, 4/5, i.e. four different 
frequencies of occurrence of stacking faults. 



shift the spacer blocks, and thus double the unit cell vol- 
ume. To modify the ratio, we added a BiS2 layer, with 
its corresponding spacer layer, at the top of 6130283. We 
therefore obtain the following stacking : 2 x BiS2 / Bi202 
/ S2 / Bi202 / 2 X BiSa / Bi202 / S2 / Bi202 / 2 x BiSs 
/ Bi202. This can be done because BiS2 layers, unlike 
S2, do not shift spacer blocks. The piling up of unit cells 
is therefore not disturbed by the added block. Similarly 
we also constructed a structure in which an additional 
BiS2 every two units cells was included. In total, we cal- 
culated the electronic structure of compounds with four 
different 82/6182 ratios : 0,1, 2/3 and 4/5. The struc- 
ture of Bi82 layers is different in Bi2082 and 6130283, 
the most noticeable difference being the angle between 
Bismuth and Sulphur atoms in the Bi8 planes, which is 
larger in the parent phase. We did calculations for both, 
and the results being extremely similar (e.g. less than 
0.1 eV difference for the height of bands), we decided to 
show the results for Bi30283's Bi82 layers, as the stack- 
ing faults appear in this compound. 

Bi20S2 and 6130283 only differ by the replacement in 
the second compound of one out of two 6182 layer by 
an 82 layer. In terms of stoichiometry, as noted by Phe- 
lan et al. [20], 612082 can be written 6106182, just as 
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FIG. 2: Band structures and density of states for the four 
simulated compounds. From top to bottom, they correspond 
to the structures from left to right in Figure 1. The two last 
compounds are simulations of stacking faults at two different 
frequencies In BI3O283. 



the parent phase of the Lanthanum compound, La06i82. 
The band structure of 6i2 082 has features close to the 
Fermi energy which are similar to the undopcd La06iS2 
phase. Of specific interest is the set of 6ismuth 6p and 
Sulphur 3p bands that are found in all the band struc- 
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tures calculated (Figure 2). Further partial density of 
states calculations enabled us to confirm that exact com- 
position. In the first case, Bi20S2, they are just above 
the Fermi level, at around 0.3 eV at their minimum. This 
is very close to the corresponding minimum in LaOBiS2 
(0.2 eV). The gap is 0.7 eV wide, twice as large as in the 
LaOBiSa case (0.37 eV). 

In 6130283, one of these bands is crossing the Fermi 
level, near the R and X points, just as in LaOo.sFo. 56182. 
At these points it interacts with two other bands. These 
are in great majority composed of S 3p states, as the 
density of states shows. They are not present in the par- 
ent phase, Bi20S2, which means they come from the 82 
layer. To confirm this assignment, we plotted the local 
density of states in real space, for the range of energies 
between -1.5 eV and eV (Figure 3). This range con- 
tains the bands we are interested in, with almost nothing 
else in. The results are quite clear : it is indeed the 
82 layer that gathers all the electron density. Moreover, 
this density is organised as a torus oriented along the z 
axis around each 8 atom. The fact that the density de- 
creases near the atom comes from the pseudopotentials, 
but we do learn that the orbitals lying at or just above 
the Fermi level correspond mainly to Px,Py electrons on 
the 82 dimers (oriented along z). The concentration of 
the electron density on the atoms rather than on the 
8-8 bond could be indicative of a weakening of the 8- 
8 interaction. Indeed, the simulation indicates there is a 
repulsive force of 5.5 eV/A for each of these atoms, along 
the 8-8 bond. 

We also plotted the local density of states in real space 
for the range of energies between eV and 1.5 eV (Figure 
3). This range contains the characteristic bands present 
in all the BiS2-based compounds. We can clearly see that 
the charges are localised in the Bi82 layer, more precisely 
in the two Bi8 planes. Indeed the 8 atoms in the Bi82 
layer but outside the Bi8 planes gather as many electrons 
in that range of energy as the 8 atoms in the 82 layer. 

The electronic structure of 6130283 with an added 
stacking fault does not differ a lot from 6130283 's (Figure 
2). However three things do change : the height of the 
Fermi level, of the 82 bands, and the 6i8 bands around 
r. The 82 bands are much closer to the lower bands, 
while the energy difference between 618 bands and the 
lower bands is preserved. A gap appears, of about 0.6 
eV, between these two sets of bands. The Fermi level is 
lower, the 82 bands not raising it as much as in 6130283. 
In consequence, the characteristic 6i8 bands around R 
and X do not cross the Fermi level. 6ut they do cross 
it around F, the bands being lower there than at R or 
X, unlike any 6182 compound encountered so far. The 
band structure of 6130283 with one stacking fault every 
two unit cell is similar. The only differences visible are 
in the density of states, well above or under the Fermi 
level. The height of the 82 and 618 bands does not seem 
to be affected. 

The key result of this study is the fact that, in spite of 
it being a nominally stoichiometric composition (no dop- 
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FIG. 3; Local density of states of 3130283 in real space for en- 
ergies in [-1.5,0] (left) and [0,1.5] (right), plotted with XCryS- 
Den [27]. Framed is a 3D plot of an isosurface centred on the 
S-S dimer. In the first case, the S2 layer gathers most of the 
density. In the second case it is the BiS planes, without the 
extra S atoms, that gather the density. Half of the atoms 
represented are not in the plane corresponding to the density 
plot. 



ing), the 82 layers in the 6130283 compound act as effec- 
tive dopants of the parent, insulating compound. They 
push a Px,Py-^i^e band of the 82 layer onto the bottom 
of the conduction band, thereby transferring part of the 
electrons of the 82 dimers to the 6182 plane. Although 
the actual band energies are not quantitatively accurate 
(the quasi-particle gap would be larger than the GGA 
band gap, on one hand, but the spin-orbit coupling would 
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push down one of the BiS2 related conduction bands, 
on the other) the quahtative fact remains, in agreement 
with experiment. Indeed, the results of our calculations 
are in agreement with the observed metallic behaviour 
of 6130283 [20]. Furthermore, the experimental obser- 
vation of an anomalously short S-S distance in the S2 
layers is also consistent with electron depletion of the S2 
layer. Indeed, a substantial residual force of 5.5 eV/A 
is observed for the S atoms in the S2 layer in the direc- 
tion of elongating the distance, which suggests a larger 
depletion than the one obtained in the calculations. 

The link between superconductivity and the crossing 
of BiS2 bands is supported by the evolution of super- 
conductivity and of band structure with stacking faults. 
The latter lower the S2 bands, causing these bands not 
to cross the Fermi level. Everything points towards this 
lowering being responsible for the loss of superconductiv- 
ity when stacking faults become numerous. 

As the band structures of the two compounds contain- 
ing stacking faults are virtually inseparable around the 
Fermi level, we can infer that their influence decreases 
slowly with its frequency of appearance. This is in ap- 
parent contradiction with the fact that superconductivity 
has been reported in a compound having around 20% of 
Bi2 0S2-like stacking, equivalent to the 82/6182 =4/5 
ratio we simulated [ ]. But when looking at the corre- 
sponding TEM data, one can see that the stacking faults 
are gathered together. The magnetization measurement 
shows the superconducting volume of the compound is 
divided by two when the Bi2082-like stacking goes from 
<10% to 20%. Therefore it is reasonable to postulate 



that the non-superconducting phase is a volume sur- 
rounding the Bi2 082-like stacking. As a consequence, 
Bi30283 is extremely sensitive to this type of disorder. 
Further refinement of purity and crystal growth com- 
bined with electronic structure measurements and mea- 
surements under pressure are required to provide further 
insight into both normal and superconducting states of 
these fascinating materials. 

Finally, the calculations of density of states in real 
space indicate that the BiS2 bands, seeming to be central 
for superconductivity in these compounds, actually come 
from the Bi8 planes. Therefore finding other structures 
having such planes could be very interesting in order to 
tune their interactions differently, and maybe transform 
their behaviour. 
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